Coherent Transport through an interacting double quantum dot: Beyond sequential 
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Various causes for negative differential conductance in transport through an interacting double 
quantum dot are investigated. Particular focus is given to the interplay between the renormaliza- 
tion of the energy levels due to the coupling to the leads and the decoherence of the states. The 
calculations are performed within a basis of many-particle eigenstates and we consider the dynamics 
given by the von Neumann-equation taking into account also processes beyond sequential tunneling. 
A systematic comparison between the levels of approximation and also with different formalisms is 
performed. It is found that the current is qualitatively well described by sequential processes as 
long as the temperature is larger than the level broadening induced by the contacts. 



I. INTRODUCTION 

The field of electronic transport through metallic and 
semi-conducting quantum dots has been a topic of in- 
tense research for over a decade. While in the early ex- 
periments the main focus was on charging effects leading 
to the phenomenon known as 'Coulomb blockade'^ and 
later on the Kondo effect in single quantum dots,- 2 in re- 
cent years the attention has shifted to more elaborate 
systems such as double quantum dot systems i 3 ' ' 5 ' 6 ' i 8 > ^ 
In addition, the study of the electronic spectrum of quan- 
tum dots (excited states) was possible as the physical size 
of quantum dots could be further reduced by improved 
lithographical methods as well as the rise of new materi- 
als such as nanotubes and semi-conducting nanowires. 10 

Already in the '90s, there were experiments^ on sin- 
gle quantum dots displaying non-monotonous current- 
voltage characteristics, for which the current decreases 
with increasing bias, leading to a negative differential 
conductance (NDC). Within the orthodox theory of se- 
quential tunneling, such effects were explained by the 
presence of excited states which were more weakly cou- 
pled to the leads than the ground states (for a given 
charge on the dot). The reason for a state dependent 
coupling could be either due to spin (Clebsch-Gordon 
coefficients) 12 or the different orbital wave- functions of 
the various statesj 13 i 14 A bias voltage dependence of the 
lead-dot tunnel coupling can lead to a weak NDC effect, 
see, e.g., Ref. [IH 

For a double quantum dot (DQD) system with the 
two quantum dots in series, very sharp current peaks 
(and corresponding NDC) was observed in experiments 
where the inter-dot coupling fl was weak. 3 Depending on 
whether two levels on different dots are 'aligned' or not, 
the current will be high (in alignment) or low (off align- 
ment). As the quantum levels in the different dots are 
shifted differently by the applied bias voltage (depending 
on the various capacitances of the system) the alignment 
condition is fulfilled at certain small ranges of the bias 



voltage, leading to the current peaks. 

In this paper we focus on quantum transport through 
a DQD. This system offers the possibility to study the 
interplay between the coherent quantum mechanical os- 
cillation inside the DQD and the influence of the coupling 
to leads. Especially, we are interested in the negative 
differential conductance caused by this interplay (rather 
than the well-known sources discussed above). In the 
regime where the inter-dot tunneling coupling dominates 
over the coupling to leads, it was found in Ref. [16| that 
NDC only occurs if the spatial symmetry of the system 
is broken, e.g. due to asymmetric coupling to leads or 
detuning of the bare level energies 4^ In another recent 
paper, B. Wunsch et al. [I?) investigated the transport in 
the opposite limit, where the inter-dot coupling is weak 
and with a small detuning of the bare level energies, i.e. 
for an asymmetric system. They found that NDC can be 
caused by a level renormalization due to the coupling to 
leads. Also in the weak inter-dot coupling regime, Djuric 
et al. [l8| found NDC even for symmetric systems for 
certain ratios between the inter-dot tunneling coupling 
and the coupling to leads. The effect was explained in 
terms of decoherence due to the coupling to leads, which 
depended on the occupation of the dot. In this article we 
show how these effects relate to each other. 

In the above mentioned works, the current was only 
calculated to lowest (first) order in the lead-dot tunnel 
coupling, so strictly speaking, the results are only valid 
in the sequential tunneling limit. However, for the issues 
addressed above, the coupling to the contacts strongly in- 
fluences not only the occupations but also affects the na- 
ture of the transport, especially its quantum-mechanical 
coherence. Thus it is not a priori clear, if the first-order 
approach is appropriate, even if the temperature is higher 
than the energy scale (line-width) due to the coupling to 
the leads. Below, we investigate the current to higher 
order in the lead-dot coupling by applying the method 
described in Ref. [l9| and compare with first order results. 
For both cases we find qualitative agreement if the tern- 
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perature is higher or comparable to the line-width due to 
coupling to the leads, but we also discuss the behavior 
for lower temperatures, where the first order approach 
becomes unreliable. 

The paper is organized as follows: in section [IT] we 
present in detail the DQD model system we consider. 
Section IIIII discusses in brief the transport models we 
are using to obtain the transport results (details are pre- 
sented in the appendix). The potential sources of NDC 
behavior are discussed in detail in section Hvl Finally we 
summarize our findings in section IVl 



II. THE MODEL SYSTEM 

We consider a double quantum dot system, where the 
spin degree of freedom has been omitted in order to sim- 
plify the analysis (in subsection llV El the double spin case 
is briefly addressed) . In a single-particle basis the Hamil- 
tonian for the system reads 



H = E a did a +E p dld + Ud\d a d^d +{ndld a + h.t 

+" E U c ll c U + {tkL d a c hL + 



k-( 



(1) 



where the first line describes the isolated quantum dot 
system with U being the Coulomb energy for occupy- 
ing both dots, fl the inter-dot tunneling coupling and 
with a. I ' j3 denoting the left /right dot. The first term 
in the second line accounts for the leads with index 
I = L/R for the left /right lead and levels counted by 
k. The last term is the lead-dot tunneling coupling. 
We parameterize the lead-dot coupling parameters tkt by 
Te(E) = 2irJ2k \tki\ 2 S(E — E kl )- Here we use the con- 
stant value Ti for \E\ < 0.95W and assume T e (E) = for 
| .El > W . For 0.95JT < \E\ < W we interpolate with an 
elliptic behavior in order to avoid discontinuities. Fur- 
thermore, we define r = Tl + Tr. The bias voltage Vt,i as 
is applied symmetrically to the electrochemical poten- 
tials of both leads, [II — — Mi? = eVbias/2, where e is the 
positive elementary charge. 

Throughout this paper we include the Coulomb in- 
teraction by considering a basis of many-particle states 
for the isolated DQD, which allows for a consistent de- 
scription of many-particle effects, see also Ref. and 
references given therein. Thus, we diagonalize the first 
line of the Hamiltonian from Eq. (fTJ) and find the eigen- 
states and the corresponding energies. Eq = is the 
energy of the empty state, E\ = p — VA 2 + 4fi 2 /2 and 
E% = p+VA 2 + 4£l 2 /2 the energies of the single occupied 
states, and Ed = E a +Ep + U the energy of the double oc- 
cupied state, where A = E a — Ep and p = (E a + Ep)/2. 
The states with energies E1/E2 are referred to as the 
bonding/anti-bonding state. 

Depending on the occupation of the dot states different 
transport regimes can be defined. Current through the 



DQD is effectively blocked if no one-particle excitation 
lies in the bias window between the Fermi levels of both 
contacts. This is known as the Coulomb blockade regime. 
Therefore, as the bias is increased, a current can flow 
through the structure whenever a one-particle excitation 
becomes energetically allowed, leading to a step feature in 
the current and a corresponding peak in the differential 
conductance. Having 4 such possible excitations (0 «-> 
1,2, and 1,2 <-> d), at most four steps can be observed 
in the TV-curve. Further steps can be seen, if spin is 
considered as welli^L 

The quantum rate equations from Refs. 22.23.24] are 
valid in the high-bias limit, i.e. if the energy difference 
between the chemical potentials in the contacts and the 
excitations exceeds both the level broadenings Ti and 
the temperature. If only one-particle states are within 
the bias window, but double occupation is forbidden, i.e., 
(Ed — E\), (Ed — E2) /!l,M-R! these equations provide 
the plateau current 



h= e r 



hW(2 + T R /T L ) + (T R /2) z + A 2 



(2) 



If, in contrast, all excitations are within the bias window, 
one obtains 



h = T- 



r L r R m 2 



/i(4r» 2 + r i r ii )(r/2) 2 + A 2 r L r iJ ' 



(3) 



These values will be compared to our calculations in the 
subsequent sections. 



III. THE VON NEUMANN APPROACH 

Our calculations are based on the von Neumann equa- 
tion for the density matrix, as described in detail in 
Ref. [H. The key idea is to use a set of many-particle 
states labeled \a),\b), . . ., with energies E a , E^ , . . . , re- 
spectively, which diagonalize the Hamiltonian of the sys- 
tem, including the many-particle interaction. (In our 
case these are the states |0), |2), \d) introduced in sec- 
tion [TTJ ) Transport occurs by tunneling of electrons with 
a quantum number k from a lead I into the system, while 
the state is changed from an iV-particle state |o) to an 
N + 1-particle state |6). The corresponding matrix ele- 
ment is Tb a (k£). In Ref. the full correlations of up to 
two particles entering and leaving the system was taken 
into account and this approach will be referred to as the 
second order von Neumann approach (2vN) in the follow- 
ing. In addition, we apply the same concept restricting 
to single electron processes, which we call the first order 
von Neumann approach (lvN). The resulting equations 
for the lvN approach are given in App. El Both the 
first and the second order approaches include the non- 
diagonal elements of the density matrix, which allows us 
to consider the regime of both weak and strong inter-dot 
coupling il. This is demonstrated in App. [Bj where our 
method is also compared with other approaches. The 
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lvN approach neglects level broadening effects (of order 
T) and is thus expected to be valid only for ksT ^>T or in 
the high-bias limit. In contrast, the 2vN approach is able 
to reproduce such effects and gives good results above 
the Kondo temperature. 19 Due to the self-consistency, 
the 2vN approach contains in addition to two-particle 
correlations also a subset of higher order correlations. In 
App.lClwe show that the result for the Anderson model is 
identical with the corresponding result obtained from the 
Real-Time diagrammatic approach in resonant tunneling 
approximation. 25 ! 26 

The lvN approach contains sums of the form 
E k TbaT;, a ,/{E k - E v + E a - i0+) (see App. De- 
composing 



E k - E v +E a - 10+ I E k - E b , +E a j 
-m5{E k - E v + E a ) , 

the imaginary part can be related to electronic transition 
rates, while the real part acts as an effective renormal- 
ization of the transition energies between different many- 
particle states. In some calculations we will neglect all 
terms resulting from these real parts (we denote this by 
"no real parts") in order to demonstrate their relevance. 

IV. SOURCES FOR NDC BEHAVIOR 

In a real experimental double-dot structure, the ap- 
plied source-drain bias Vbi as does not only determine the 
electrochemical potentials in the leads, but it will also 
shift the dot level energies by polarization. The amount 
of these shifts depends on the details of the various dot 
capacitances and can be taken into account by lever arm 
factors X a , \p for the respective dot levels-. In addition, 

if gates are present, the gate voltages V^"{f can also shift 
the respective dot levels with efficiency factors r\ a ^r\p. 
Therefore, the voltage dependence of the dot level ener- 
gies can be written as 

E a = E° a + \ a ^p-r, a eV g l tc , (4) 

Ep = ^-Xp^-npeV^ (5) 

with Ef being the equilibrium level of the energies. This 
allows for an independent control of Vbias, the level dif- 
ference (detuning) A = E a — Ep, and the average level 
p = (E a +Ea)/2. In the following we set p = 0, meaning 
that the dot states are at equal energetic distance from 
the equilibrium Fermi level. In Fig. [T] we show the cur- 
rent calculated with the lvN approach as a function of 
Vbias and A for different inter-dot coupling strengths SI. 

In a real experiment, the current voltage characteris- 
tic corresponds typically to a line in the (Vbias, A)-plane. 
With zero gate voltages, one has 

A = E° a -E$ + (X a + Xp)^, 



(a) fi = r/10, k B T = T 




eV him [F] 



(b) n = r/4, k B r = r 
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(c) a = 2r, k B T = r 
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FIG. 1: Current versus bias voltage and detuning A = E a — 
E/3 for (a) n = T/10, (b) Q = r/4, and (c) SI = 2F using 
the first order approach (lvN) including the real parts. The 
other parameters are E a + Ep — 0, Tl = Fr = T/2, ksT = T, 
U = 10r and W = 35F 

so one follows a straight line with positive slope, see, e.g., 
the dashed line in Fig. [IJa). For a sufficiently large slope 
(X a +Xp)/2, we observe first an increase of current for low 
bias, and then a decrease of the current as the levels move 
out of resonance with increasing A. This is the standard 
NDC effect induced by electrostatic polarization. 

In addition to the above 'trivial' effect we can identify 
two further scenarios for NDC for a fixed detuning A4i, 
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FIG. 2: Current versus bias voltage for different detunings 
A. All other values like in Fig. [T] The dashed lines are the 
lvN results including the real parts, the dashed-dotted lines 
are the lvN results without the real parts (noR), and the full 
lines are the 2vN results. The values from Eqs. Q2I3J I for A = 
and A = ±r/4 or A = ±2r (which do not depend on the sign 
of A) are shown on the ?/-axis. We only show positive bias 
here, as the negative bias result corresponds to the results 
with the opposite sign of A for the symmetric coupling to 
contacts considered here. 



which will be discussed in detail below. Firstly, we no- 
tice that for small and intermediate £1 [Figs. [Ha,b)] the 
current peak (red region) is shifted to negative A with 
increasing bias < eVbi as < 20r = 2U. Thus, the cur- 
rent drops with Vbias if A > is kept constant, see also 
the dashed lines in Fig. HJa,b). A more detailed analysis 
is given in Sec. IIV Bl 



Secondly, the height of the current peak (at fixed bias) 
drops when the bias voltage exceeds 2U = 20T in Fig[TJa) 
where the inter-dot tunneling coupling is much weaker 
than the coupling to leads. This provides NDC around 
A « as discussed in detail in Sec. IIV Dl With increasing 
inter-dot coupling this effect vanishes as seen in Figs.[]Jb- 
c). 



A. First order versus second order 

In Fig. [5] we provide a systematic comparison between 
the lvN approach (dashed line) and the 2vN approach 
(full line) for different values of the detuning A, corre- 
sponding to cuts along horizontal lines in Fig. [TJ We find 
that both approaches are in good qualitative agreement 
both for large and small values of the interdot coupling 
at the considered temperature. This shows that the lvN 
approach works well even for the moderate temperature 
fcflT = Tl +r^. We observe small discrepancies close 
to current steps, where the broadening is underestimated 
due to the neglect of line-width broadening in the lvN 
approach. As expected, these discrepancies are strongly 
enhanced if the temperature drops below the level broad- 
ening as shown in Fig. [^b) for fc^T = (Tl + Tn)/5. 



B. NDC due to level renormalization 

Let us now focus towards the bias range eVbias < 2C7, 
where the double occupied state does not yet contribute 
to the current. In this regime we observe a significant 
shift of the current peak from its naively expected posi- 
tion at A = for small and intermediate il [Figs.[TJa,b)]. 

In Ref. [I?] Wunsch et al. considered a DQD includ- 
ing spin. They calculated the transport using a first- 
order diagrammatic real-time transport approach (see 
e.g. Refs. [HH, restricting themselves to the limit of 
small inter-dot coupling il <C T, where electronic states 
localized on the single dots form an appropriate basis. 
For positive A they observe pronounced NDC similar to 
Fig. [2Ja,b). This was explained in the following way: 
The energy levels E a , Ep are renormalized due to the 
couplings to the contacts. For finite U this renormaliza- 
tion is strongest if the levels are close to the chemical 
potential. Now, the localized states used in Ref. Il7| cou- 
ple mostly to the nearest lead and thus the renormaliza- 
tion differs for both levels at finite bias. This provides 
a bias-dependent renormalized A c g and the maximum 
of the current occurs at A e g = rather than A = 0. 
This effect docs not occur for U — because in this case 
the renormalization does not depend on the location of 
the chemical potential (this becomes obvious in a Green 
function treatment providing the exact result for U — 0) . 

Our results in Figs. [l^a,b) are in full agreement with 
these findings. In particular, the shift of the current peak 
position can be directly attributed to a bias-dependent 
renormalization A — > A e gf, where A e g > A for positive 
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(a) n = T/10, k B T = T/5 




C. Symmetry with respect to the sign of A 
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FIG. 3: Results for low temperature k B T — F/5 and H = 
r/10. All other parameters as Figs. 11121 Upper panel: Results 
of the lvN approach. Lower panel: Comparison of lvN, 2vN 
and the lvN without real parts (noR) approaches. The lvN 
strongly exaggerates the features at the current steps for this 
low temperature in comparison to the 2vN approach. The 
lvN without real parts is qualitatively incorrect. 



bias. This renormalization is also reflected in the magni- 
tude of the current: For eVbi as < 2U, Eq. ([2]) suggests the 
plateau current I\. In contrast, no such plateau is seen 
in Fig.[2ja,b) for small and intermediate tunnel coupling. 
Furthermore, the calculated currents are smaller than the 
respective value of Ii (|A|) for A = and A = r/4, while 
they exceed 7i(|A|) for A = -r/4. This can be fully 
attributed to a bias-dependent renormalized A e ff > A 
which should be used in Eq. @. 

Fig. [He) shows that the effect vanishes for strong inter- 
dot coupling fl > r, a regime where the approach of 
Ref. [13 fails. In addition, the quantitative agreement 
between the lvN and 2vN approach shows that the shift 
of the resonance condition is not an artefact of a first 
order tunneling approach (such as lvN or the approach 
used in Ref. Il7l ) but persist even if higher order tunneling 
processes are taken into account, which, e.g., cause line- 
width broadening. 



The quantum rate equation results Eqs. (|2I3|) only de- 
pend on the absolute value |A|. In contrast, Figs. ll|2l ex- 
hibit a strong asymmetry with respect to the sign of A. 
For the cases f2 < T (Fig. [5] (a,b) ) the reason lies mostly 
in the level renormalization discussed above. This can be 
deduced from the fact that the curves without real parts 
(noR) that neglect the level renormalization are actually 
approximately symmetric. However, for CI > T [Fig. [5] 
(c)] all numerical approaches give an additional step at 
eVbias/2 ~ Ei + U for positive A = 2r, which is absent 
for negative A. 

For the case of large ft the "molecular" bonding and 
anti-bonding states provide an appropriate description of 
the single-electron states of the DQD (see also the com- 
parison in App. lB]) . Because of the large (positive) detun- 
ing A > r, ksT the wave function of the bonding state 
has much more weight on the right dot (interface to the 
collector lead) than on the left dot. Correspondlingly, 
the anti-bonding state has more weight on the left dot 
(interface to the emitter electrode). This asymmetry in 
the wave functions leads to an effective asymmetry of 
coupling of these states to the leads, e.g. for positive A, 
the anti-bonding state is strongly coupled to the emit- 
ter, whereas the bonding state is strongly coupled to the 
collector lead. 

The eigenenergies E\$ of the bonding and anti- 
bonding state are given by =F"\/A 2 + 4f2 2 /2 = ^fy/5T (see 
section Hi]) . As E\ = — y5r is the only state with neg- 
ative energy, it is the overall ground state of the DQD. 
This means that in equilibrium (zero bias) the bonding 
state is mainly occupied, whereas all other states have 
only small occupation probability. The first current step 
appears at a bias eVbi as /2 ~ \E\\ ~ 2.2T when the DQD 
can be emptied. At the same bias the anti-bonding state 
can be occupied by tunneling of an electron from the 
emitter into the empty DQD. Because of the coupling 
asymmetry, above a bias eVbias/2 ~ \E\ \ there is a large 
change in the average occupation of the DQD states: the 
bonding state is strongly depleted due to the good cou- 
pling to the collector lead, whereas the anti-bonding state 
is now favorably filled by electrons tunneling into the 
DQD from the emitter lead. Thus, on the first plateau 
the DQD is most of the time in the anti-bonding state. 

As the bias is further increased, the doubly occupied 
state comes within energetic range when the bonding 
state can also be filled in addition to the anti-bonding 
state. This happens at eVbi as /2 = E\ + U = (— \/5 + 
U)T w 7.8r, which is where the second current step 
sets in. Now, the anti-bonding state loses occupation 
in favor of the double occupied state and the bonding 
state. Finally, at even larger bias eVbias/2 = E% + U = 
(V5 + 10)r ps 12. 2r, the doubly occupied state can be 
populated by an electron tunneling in from the emitter, 
even if the DQD is previously in the bonding state. This 
leads to the (weak) third current step for the curves cor- 
responding of positive values of bias and A. 
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In contrast, for negative A and positive bias, because 
of the reversed spatial asymmetry of bonding and anti- 
bonding wave functions, the DQD remains mostly in 
the bonding state, even though the bias has exceeded 
~ | -Ei | ~ 2.2F. The occupation of the anti- 
bonding state remains very small in the range 2.2T < 
eVbias/2 < 12. 2r. Therefore, the middle current plateau 
is strongly suppressed and smeared out by the thermal 
(and linewidth) broadening. 

The quantum rate equation results Eqs. (|2l3p corre- 
spond to the current values of the first and third cur- 
rent plateau. To capture the middle plateau one would 
have to account for the possible transitions from the anti- 
bonding state to the doubly occupied state in the relevant 
bias range. This was not done in Ref. ITsi consequently 
the middle plateau is never observed in their Fig. 6 even 
for the case f2 = V. 



D. NDC due to decoherence 

Now we focus on the behavior around eVbias ~ 2U, 
where the double occupied state enters the window be- 
tween the left and the right chemical potential. Fig. Ufa) 
shows that the current peak is significantly larger for 
eVbias < 2Z7 (region 1) than for eVbias > 2{7 (region 2) for 
small il. Just the opposite holds for large O, see Fig.QJc). 

This drop of current has been addressed by Djuric et 
al*£, using the quantum rate equation formalism devel- 
oped in Ref. [2?]. They restricted the analysis to a fully 
symmetric system, i.e. A = 0, and did not include any 
real parts. In this case Eqs. (12l.'>l) provide I\ > I 2 for 
H < v / FlF^/2, and h < h for O > VTlFr/2. Nu- 
merically, they observe a smooth interpolation between 
these plateau values upon variation of Vbias, similar to 
the result of the lvN approach without real parts (dot- 
dashed line) in Fig. [5] This gives rise to NDC around 
eVbias ~ 2t/ for weak inter-dot coupling ft < ^/TlTr/2. 

While the observed increase of current for £1 > 
v / TlT t r/2 can be easily attributed to the opening of a 
new current channel, the NDC for Q, < v / TlT T r/2 is less 
straightforward. In the limit of small il, Eqs. (|2I3[) read 

h « ^ 2 L(A,r fl ), I 2 ^^-Q 2 L(A,T R + T L ) 
h h 

with the Lorentzian L(A, 7) = A2+ ^ y / 2 ) :i • This is just the 
expression from Fermi's golden rule for sequential tunnel- 
ing between the localized dot states, which is limiting the 
current for weak coupling, see Ref. l28l . The broadening 7 
is given by the dephasing of the coherence between the lo- 
calized dot states^. For the second plateau I2 the inter- 
dot Coulomb repulsion does not play a role and we get 
the broadening 7 = Tr + T^. However, for eVbias < 2C/, 
the left contact cannot add an electron to the system if 
one electron is already present in the DQD. Thus, the 
dephasing of the coherent transitions between the dots 
is only due to 7 = Tr. Therefore, in Fig. Q^a), tne cul '~ 
rent peak (at fixed bias voltage) is higher and narrower 



(as a function of A) for eVbias < 2Z7 (first plateau) in 
comparison to eVbias > 2C7 (second plateau). 

As displayed by Fig. [2] the transition between the 
plateau values for A = is much more complex than 
the smooth transition suggested by Ref. |l8|. Both the 
lvN and the 2vN approach do not reach a stable value 
at the first plateau even for U being by far the largest 
energy. Instead the current drops over the full length 
of the platea u, g iving a much weaker NDC than pre- 
dicted by Ref. Uj£ The plateau value I2 is indeed reached 
in region 2 for Vbias — * 00, but due to the level renor- 
malization effect the transition is much broader than if 
it was only given by temperature (or even a combina- 
tion of temperature and line- width broadening T). For 
^ = \/TiXr/2 (Fig. [2^b)) the current does not reach the 
plateau value I\, and at the transition between the two 
regimes (eVbias — 2J7 = 20r) a dip in the current is ob- 
served. For Q, = 2T (see Fig. [2c)) where I\ < 1% and no 
NDC should occur due to decoherence, we observe that 
the plateau values are indeed reached in both regimes. 
While Fig. [2fa) shows very non-monotonous behavior, 
Fig. [U shows that the current varies more continously in 
the (A, Vbias) plane. This shows that the main difference 
between our lvN approach and the result of Ref. [H is 
the level renormalization discussed above resulting in an 
effective A ff , which is not contained in the quantum rate 
equation formalism of Ref. I27I . 

So far we have only considered a fixed value of the 
Coulomb repulsion U = 10T. Fig. Q] shows results for 
different values of U at weak coupling £1 = T/10 and 
zero detuning A = 0. It should be noted that even for 
U being by far the dominant energy (e.g. U = lbT) 
no clear plateau value is observed, but instead a slow 
cross-over between the two regimes. Again, relatively 
small differences between the lvN and 2vN approach 
are observed, though in general the 2vN approach re- 
duces somewhat the level renormalization effect and leads 
to more plateau-like current-bias characteristic than the 
lvN approach. Here, a particular surprising feature is 
the fact that a small (unphysical) NDC remains even for 
U = in the lvN approach (see also the analytical result 
in App.[B]), which is however absent in the 2vN approach 
(in full agreement with the transmission result by Green 
functions). 

We conclude that the NDC-scenario due to decoher- 
ence outlined in Ref. [H is strongly modified by the level 
renormalizations addressed in section IIVB1 Also note 
that the high-bias limit I2 is only reached for very high 
bias in the case of weak or intermediate inter-dot cou- 
pling f2 as shown in Fig. [2] 



E. Transport with both spins 

While we restricted ourselves to the case of spinless 
fermions before, we now take into account both spin di- 
rections in both dots, thereby accounting for sixteen dif- 
ferent many-particle states. We add an intradot Coulomb 
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FIG. 4: Results for A = 0, fi = T/10, and different values of 
U. All other parameters as in Figs. 1 1 121 



n = r/io, k B T = r 




10 20 30 40 50 



F bias [r] 

FIG. 5: The current taking into account both spin directions 
and an intradot Coulomb repulsion of (/intra = 30F calculated 
by the lvN approach. All other parameters as Fig. []Ja) . 

repulsion of (/intra = 30r between the different spin states 
within each dot, while the interdot Coulomb repulsion 
U = I0r is assumed to be independent on the spin di- 
rection. The result from the lvN approach is shown in 
Fig. which is qualitatively similar to Fig. [Ha) where 
spinless electrons were considered. However, quantita- 
tively, the features of level renormalization and NDC due 
to coherence are enhanced. 



V. DISCUSSION AND SUMMARY 

We have studied transport trough a double quantum 
dot system and investigated various sources of negative 
differential conductance (NDC) potentially observable in 
experiment. Our method reproduces the basic features of 
Refs . 11711181 . which constitute certain limits of our full nu- 
merical approach. In particular we can treat all values of 



inter-dot coupling and detuning of dot levels within the 
lvN approach. Further effects of line-width broadening 
due to higher order tunneling events can be taken into ac- 
count by the 2vN approach. The NDC due to level renor- 
malization (Sec. IIVB[) as introduced by Ref. [l7J using 
a strictly sequential tunneling scenario is only quantita- 
tively modified by the higher order tunneling events (even 
at temperatures comparable to the line-width broaden- 
ing), thus showing a surprising robustness. The NDC 
due to decoherence^ (Sec. HVDj) is clearly seen in the 
heights of the current peaks in the bias-detuning plane, 
see Fig. Q] Nevertheless this effect is strongly masked by 
the level renormalization effects if a constant detuning A 
is considered. Both effects show the relevance of a consis- 
tent treatment of first-order tunneling terms, which can 
be achieved by the lvN approach discussed here. Com- 
parison with the higher-order 2vN approach^ for a large 
variety of parameters shows that the validity of these 
first-order approaches is only restricted to the tempera- 
ture being larger than the level broadening. 



Acknowledgments 

This work was supported by the Swedish Research 
Council (VR) and the Villum Kann Rasmussen fond. 



APPENDIX A: FIRST ORDER VON NEUMANN 
APPROACH 

We denote the many particle states by |o), \b) . . . with 
energies E a ,Ef,,..., respectively, and we use the conven- 
tion, that the particle number follows the position of the 
letter in the alphabet, i.e. \a) is an iV-particle state and 
| b) is an N+ 1-particle state. The tunnel matrix elements 
Tb a (ka£) for a transition from the state |o) to \b) by enter- 
ing of an electron from the contact I € {L, R} with spin 
a and momentum k can be directly related to the single 
particle tunneling Hamiltonian as given in Eq. JT]), see 
App. A of Ref. |l9|. The key quantities are the elements 
of the reduced density matrix Wb'b = Trace {\b) (b'\p}, 
where the diagonal elements are the probabilities to find 
the respective many-particle state and the off-diagonal 
elements refer to correlations between the many-body 
states induced via coherent tunneling processes to and 
from the leads. 

The equation of motion for Wb>b is derived from the von 
Neumann equation for the density matrix, see Eq. (11) 
of Ref. fl9L and depends on the current amplitudes 

4> c b{kai) — Trace||&)c^(c|p|. These 4> c b(ka£) are 

themselves determined by an equation of motion, see 
Eq. (7) of Ref. [H In the first-order approach, all terms 
containing both k and k' (correlations between two tran- 
sitions) are neglected and Eq. (10) of Ref. [H has the 
solution 
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= ~f dt'e^+^-^+^-t'y 1 (^2T cb ,{k)w Vb {t')h -Y^ Wcc ,(t')T c , b (k)(l - fk)j (Al) 

Now we neglect the time dependence of Wybit') m the kernel of the integral (the Markov limit) and set Wybit') = 
Wb>b(t), which allows us to perform the integral. Inserting this into the equation of motion for %j, we obtain 



ih^-Wbb' ={E b - E b >)w b b> 



£«' ™ aa >TZ, al (ka£)h(E k ) - £ b „ T*„ a (ka£)wb^(l - fe(E k )) 



E k - E v + E a - i0^ 

a.kat 



Eg' Tba>(kae)w a , a f e (E k ) - J2 b » w bb »T b „ a (kae)(l - fifth)) , , , 



,k<j£ 



E k - E b + E a + 10+ bay ' (A2) 



, ST rr*ri, ^J2b" T cb"(kcr£)wb" b 'fe{E k )-J2c' w cc'T C ' b '{kae)(l-fe(E k )) 

+ W^ kal) E k -E c + E b , + m 



X- E b » w w ,T c * b „(ka£)ME k ) - E c , T* b (ka£)w c , c (l - fifth)) ^ 
~ \ f E k -E c + E b -m Tcb ' ^ ' 

The current from the left lead into the system is given by 

J L = -h< V -= :/ I > T^Lwrn!^)- > «v,,r,,,OH(1 - A (£,,)) 1 J- . (A:-!) 



£ Ek Ji c { +E b + iO+ [^T cV {kaL)w b .bf L ft k )-^w cc ,T clb {kam 



If we restrict to diagonal elements P b = Wbb, these equa- 
tions reduce to the standard master equation 30 ^ formu- 
lated in a many particle basi s 31 ' 32 ! 33 . 

The Redfield kinetics 3 ^ has been recently used to 
derive a similar set of equations 3 ^. We can re- 
cover these equations, if we approximate w b > b (t') — 
w b 'b{t)e~^ E "'~ Eb '> (t '^ /h in Eq. fAT]) . In this case the 
integrations can be performed as well, but we obtain 
slightly different denominators for the nondiagonal ele- 
ments in Eq. (|A2[) . While the approximation w b > b (t') = 
Wb'b(t) becomes exact in the stationary state, which we 
consider here, the behavior Wb>b(t') oc Q- 1 ( E b<- E b)' t l h [ s 
suggested by the linear term in the equation of motion.— 
At the moment, we have no direct indication, which con- 
cept is more appropriate. However, we did only find mi- 
nor numerical differences and the qualitative features are 
identical for all issues discussed in this article. In par- 
ticular, both approaches can yield negative probabilities 
Wb b , a well-known problem of Redfield kinetics 3 ^. 



APPENDIX B: COMPARISON OF DIFFERENT 
APPROACHES IN THE NONINTERACTING 
LIMIT 

Now we consider the double dot model without spin in 
the noninteracting limit U — 0, and set E a = Ep = 0, 
= r# = r/2 as well as W — > oc. In this case we can 



solve most approaches analytically, allowing for a better 
understanding of the structure of the various approaches. 

1. Transmission formalism and 2vN approach 

As a bench mark, for the noninteracting case the cur- 
rent can be evaluated exactly via the transmission for- 
malism (see e.g. Ref. |37| ) 

I = 1 knj dET ( E )iM E ) - ( B1 ) 

with 

r 2 ^ 2 

T ^ = 4[{e - n) 2 + r 2 /i6][(s + n) 2 + r 2 /i6] ^ B2 ^ 

where the wide band limit is applied. We obtained nu- 
merically the same result from the 2vN approach for all 
parameters checked. 

If r <C ksT, we may replace the peaks in the trans- 
mission function by (5-functions which provides us with 




The prefactor for fl <C T is chosen such, that the integral 
over E agrees with the full transmission function for all 

Q,r. 
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2. Master equation 

The master equation in the many-particle state o 31 ! 32 
can be derived by setting ww = Pb&w m the lvN ap- 
proach, resulting in the current 

/master = ^ [h (A) " /fi(O) + /i(-fi) - /*(-«)] 

8/1 

(B4) 

which exactly equals the bench mark (|Bip in the limit 
r < fc s T,ft, see Eq. (jB3)l . 



3. Quantum rate equation 

Going beyond the master equation, correlations be- 
tween different states can be taken into account. The 
dephasing of these correlations is frequently treated in a 
Lindblad form 3 ^, see, e.g., Refs. [HJ^Hil]^. They can 
be derived in different ways and the name "quantum rate 
equation" is frequently used. Formulated in a basis of the 
localized states, see e.g. Ref. or Eq. (36) of Ref. [27], 
we find the result 



O 2 



^quantum rate 



■(/l(O)-MO)), (B5) 



which matches perfectly the bench mark result in the 
limit £1 < T < k B T, see Eq. (lB3|) . In addition, it is cor- 
rect in the high-bias limit, (— [i R ) 3> as proven 
by Gurvitz and Prager.— 



4. lvN approach 



Finally the lvN approach provides 

t - 1 r " 2 [h(n)-f R (n) + f L (-n)~f R (-n)} 

lvN ~HT 2 /4 + An 2 2 

a nr 2 



(B6) 



where 



A= l -jyE[f L{ E)-f R{ E)]v{-^ 



n E-n 



is a small contribution from the real parts. If the real 
parts are neglected (A — 0), the result matches the bench 
mark result as long as k B T > T. This indicates that 
the inclusion of the real parts may not be appropriate 
in the noninteracting limit for a first order approach in 
the coupling T, see also the little peak at eVbias = 2r 



for U — in Fig. [4] However, the comparison with the 
2vN approach (which does not display the spurious peak) 
indicates that the real parts cover the essential physics 
in the interacting case. It is interesting to note, that 
the Redfield kinetics provides exactly the same analytical 
result (|B6p . so that both approaches exhibit the same 
problem in the noninteracting case. 

APPENDIX C: COMPARISON WITH 
REAL-TIME DIAGRAMMATIC APPROACH 

The 2vN approach allows for an analytic solution for 
the single dot model with spin (the Anderson model) 
with infinite Coulomb repulsion. We consider a spin- 
degenerate level with E^ = E± — Ed- Analogously to 
section III of Ref. [l9|, we define 

k 

B a , {E) = B L a . fl (E) + B%(E) 
T e (E) = 2nJ2S(E - E k )T 2 {k) 

k 

r(E) = r L (E) + r R (E) 

where we assume that the couplings to the contacts 
T T0 (fc t I) = T l0 {k 11) = T e (k) do not depend on spin. 
Then we obtain 



iH-Bl fi (E) = (E d -E + £(25) + E'(J5))i£ (25) 



2tt 



Fe(E) f dE , B^ (E') + fj(E)Bl (E') 



2ir J E-E' + i0+ 

where a denotes the spin opposite to a and 

k f e E k - E k , + i0+ 



S(^)=E 



^)=£^' (fc ' )2 



f E k - E k , + i0+ 



k'f 



as well as 



2Z^J dEB hQ (E)+B L0 (E)} (C3) 



at 



With the Ansatz B^. = B^.q = B(E) and -j = w\.\ 
we have the stationary solution 



10 



w ,o[T L (E)f L (E) + T R (E)f R (E)} w^l{E)(1 - f L (E)) + T R (E)(l - f R (E))} 

As 

3{£(£) + &(E)} = -[T L (E)(1 + f L (Ej) + T R (E)(1 + f R (E))}/2 (C4) 



we find that there is a solution B(E) which is purely real (like in the spinless level case). Inserting into Eq. (|Clj) gives 
the stationary state 



r, „ T , x rJl + f L (E)}(E d - E + Y,(E) + Y,f(E))B(E) 

(Ed — e + m + e^))^) = ^r^f^ + r^ + U)) 

r£r fl (/ fl (£0 - f L (E)) 



2ir[T L (E)(l + f L (E)) + T R (E)(1 + f R (E))} 



(C5) 



where we used wo ; o + 2-uv<j = 1. With Eq. (|C4j) we have 
.fpL rnl _ r L (E)T R (E)(f R (E)~f L (E)) 

Finally the particle current (including both spin direc- 
tions) is given by: 

T = i / r L (£;)r fl (g)(/ L (£;)-/ fl (£?)) 
L ^7 ^i^-s + e^ + e/^)! 2 



This result fully agrees with the result from real-time 
perturbation theory in the so-called resonant tunneling 
approximation (containing a resummation of diagrams 
beyond second order perturbation theory), see Eq. (4.61) 
of Ref. 39. This indicates that the 2vN approach contains 
an equivalent set of higher than second order tunneling 
processes due to the self-consistency in the equation of 
motion for </> c fc(fc). 
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The results in Ref. [l6j were derived for a double dot in- 
cluding spin, but the conclusion can be transferred to the 
spin-less case. 

A constant A can be realized by minimizing the polariza- 
tion of the dot levels from the leads (vanishing A) or by 
compensating with appropriate gate voltages proportional 

tO Vbias- 

For the 2vN-approach, the same ambiguity in closing the 
equations arises. In this case, the analytical result for the 
single level problem only matches the Green's function re- 
sult, if the time dependence is neglected, which is the ap- 
proximation scheme used in Ref. IiqI . 



